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ABSTRACT 

We present timing analysis for a Rossi X-ray liming Explorer observation of Cygnus X- 1 in its hard/low state. 
This was the first RXTE observation of Cyg X-l taken after it transited back to this state from its soft/high state. 
RXTE’s large effective area, superior timing capabilities, and ability to obtain long, uninterrupted observations 
have allowed us to obtain measurements of the power spectral density (PSD), coherence function, and Fourier 
time lags to a decade lower in frequency and half a decade higher in frequency than typically was achieved with 
previous instruments. Notable aspects of our observations include a weak 0.005 Hz feature in the PSD coincident 
with a coherence recovery; a ‘hardening’ of the high-frequency PSD with increasing energy; a broad frequency 
range measurement of the coherence function, revealing rollovers from unity coherence at both low and high 
frequency; and an accurate determination of the Fourier time lags over two and a half decades in frequency. 
As has been noted in previous similar observations, the time delay is approximately proportional to and 
at a fixed Fourier frequency the time delay of the hard X-rays compared to the softest energy channel tends 
to increase logarithmically with energy. Curiously, the 0.0 1-0.2 Hz coherence between the highest and lowest 
energy bands is actually slightly greater than the coherence between the second highest and lowest energy bands. 
We carefully describe all of the analysis techniques used in this paper, and we make comparisons of the data 
to general theoretical expectations. In a companion paper, we make specific comparisons to a Compton corona 
model that we have successfully used to describe the energy spectral data from this observation. 

Subject headings : accretion — black hole physics — X-rays: binaries 


1. introduction 
1.1. Plan of Paper 

Recently, we have presented the results of a spectral anal- 
ysis of a 20 ks Rossi X-ray Timing Explorer (RXTE) obser- 
vation of perhaps the most famous of the galactic black hole 
candidates (GBHC), Cygnus X-l (Dove et al. 1998, hereafter 
paper I). The results of this analysis are summarized in §1.3 
below. Here in this paper we focus on the timing analysis of 
our RXTE observation. Specifically, we present the Cyg X-l 
power spectral density (PSD), the coherence function, and the 
time delays between soft and hard lightcurve variability. We 
discuss in general terms the analysis techniques that we have 
used. As the cross-spectral techniques utilized for determining 
time lags and coherence are unfamiliar to many researchers, we 
shall adopt a somewhat “pedagogical tone” at times. Practi- 
tioners in the field themselves have not always properly defined 
error bars and noise limits (cf. §3.1, §4.2). Cross spectral tech- 
niques are used increasingly in the study of X-ray binaries, and 
they have benefits relative to the more conventional cross corre- 
lation function, particularly in measurements with large signal- 
to-noise ( S/N ). We hope that the discussion here is of general 
use. 

We begin this paper with a review of the parameters of the 
Cygnus X-l/HDE 226868 system in general (section 1.2), and 
the modeling of the Cyg X-l high energy spectrum in specific 
(section 1.3). The rest of the paper is devoted to studying the 
variability of our RXTE Cyg X-l data. In section 2 we define 
the periodogram for one lightcurve, i.e. the PSD, and describe 


how to account for instrumental effects. We present our re- 
sults for Cyg X-l in five different energy bands. Comparing 
the lightcurves measured simultaneously in these different en- 
ergy bands can reveal further information on the nature of the 
underlying physical processes. Section 3 describes our compu- 
tations of the coherence function, while section 4 describes our 
measurements of the time delay between the variability compo- 
nents of different energy bands. Coherence measures the degree 
of linear correlation between time series as a function of time 
scale, or equivalently of Fourier frequency. In both of these 
sections, we illustrate the techniques, estimates of noise lev- 
els, etc., with our RXTE data for Cyg X-l . More importantly, 
we apply the techniques to the Cyg X-l data in order to help 
constrain physical models of this system. In section 5 we inter- 
pret our Cyg X-l results in terms of simple, physical toy mod- 
els. Specific results and comparisons for the Comptonization 
model of paper I are presented in a companion paper (Nowak 
et al. 1998, hereafter paper HI). 

1.2. Parameters of the Cygnus X-l/HDE 226868 System 

Cygnus X-l is one of the most firmly established persistent 
GBHC. Discovered in a rocket flight in 1964 (Bowyer et al. 
1965), this object was one of the first X-ray sources known, 
although its exact nature was not clear until its optical identi- 
fication with the O-star HDE 226868 (Hjellming 1973). Since 
then, Cyg X-l has been the subject of observations by almost 
all X-ray astronomy missions. In-depth reviews have been pub- 
lished by Oda (1977) and Liang & Nolan (1984). Table 1 lists 
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the most important parameters of the system. 

The black hole mass given in Table 1 is lower than the usu- 
ally quoted value of 16±5Af© determined by Gies & Bolton 
(1986). These authors determined the mass by using a relation- 
ship between the rotational broadening of the absorption lines 
of HDE 226868 and the radial velocity, and by assuming that 
these two values depend only on the mass ratio. By modeling 
the light curve they found Af = 33 ±9 A/© for HDE 226868. The 
black hole mass of 16 Af© was then determined using the mass 
function given in Table 1 . This fairly high value has propagated 
into the literature, and is quoted in most reviews on black holes 
(e.g., Tanaka & Lewin 1995). It should be noted, however, that 
the method used by Gies & Bolton (1986) is nonstandard, and 
that their mass estimate for HDE 226868 disagrees with those 
obtained by other methods. Masses obtained from photometric 
modulations are in the range of 1 6 ± 2 A/© for HDE 226868 and 
8 ± 3 Af©for the black hole (Balog, Goncharskil & Cherepash- 
chuk 1981; Hutchings 1978). Masses following from a determi- 
nation of the log# value of the O-star atmosphere typically re- 
sult in masses of around 15 A/© for HDE 226868 (Herrero et al. 


For our purposes, the main parameters of interest in Table 1 
are the compact object mass, system luminosity, and system 
distance. However, few reviews have included all the poten- 
tial parameters of interest. Furthermore, many reviews quote 
secondary sources. For these reasons, we present Table 1. 

1.3. Spectral Results 

The X-ray spectrum of Cyg X-l and other galactic black 
hole candidates in the hard state can be represented as the sum 
of a soft thermal component and a hard, non-thermal, power 
law component with a high energy cutoff at about 150keV. 
The spectrum is further modified by reprocessing features, in 
particular a weak iron fluorescence line at 6,4 keV and possi- 
bly a Compton reflection hump at energies around 30keV. All 
four components are usually explained in the context of accre- 
tion disk corona models (Hua & Titarchuk 1995; Svensson & 
Zdziarski 1994; Coppi 1992; Sunyaev & Titarchuk 1980; Sun- 
yaev & Triimper 1979). These models assume that the thermal 
radiation originates in an accretion disk and is Comptonized 
in a hot plasma, the “accretion disk corona” (ADC). Comp- 
tonization results in a non-thermal power-law spectrum. The 
iron line and the reflection hump suggest that part of the hard 
radiation from the ADC is reprocessed by the comparatively 
cold (150eV) accretion disk. 

The exact geometrical configuration of the accretion disk 
and the corona is unknown, and is one of the major subjects 
of black hole X-ray astronomical research. Due to the simi- 
larity between the hard state X-ray spectrum of galactic black 
hole candidates and active galactic nuclei (AGN), most previ- 
ous work has assumed a “sandwich configuration” where the 
corona is situated above and below the accretion disk (Pouta- 
nen, Svensson & Stem 1997; Hua & Titarchuk 1995; Haardt 
& Maraschi 1993 and references therein). The evidence for a 
sandwich configuration in GBHC, however, is not as strong as 
in the case of AGN. In previous work, some of us have pointed 
out that the efficiency of Compton cooling might prevent the 
sandwich corona from getting hot enough to generate the hard 
X-ray spectrum observed in the hard state of Cyg X-l (Dove, 
Wilms & Begelman 1997; Wilms et al. 1997; Gierliliski et al. 
1997; Poutanen, Krolik & Ryde 1997). Due to the high geo- 


1995; Sokolov 1987; Aab et al. 1984, and references therein). 
Since the work by Herrero et al. (1995) represents the “state of 
the art” in the theory of hot stellar atmospheres including wind 
loss, we present their value for the mass of HDE 226868. We 
consequently use a smaller black hole mass in this paper. This 
mass is still high enough, however, to point without doubt to 
the presence of a black hole in the system. 

A problem intimately connected with estimating the black 
hole mass is the determination of the orbital inclination. The 
value of 35° in Table 1 is presented without a statement of 
its uncertainty. Inclination measurements are intrinsically im- 
precise because the inclination has to be determined from dif- 
ficult optical polarization measurements and consequently the 
values found for the inclination are scattered ffom 26° to 67° 
(Bochkarev et al. 1986; Ninkov, Walker & Yang 1987a; Ninkov, 
Walker & Yang 1987b; Dolan 1992, and references therein) 
with most authors deriving rather lower inclinations. For these 
reasons we use i = 35°, the value adopted by Herrero et al. 
(1995) in their mass determination shown in Table 1 . 


metrical covering factor of the corona, about 50 per cent of all 
hard radiation from the corona enters the cold accretion disk, 
where it is either reflected or thermalized. The soft photons 
produced this way again enter the corona where they efficiently 
Compton cool the plasma. For a given Thomson optical depth 
of the corona, there thus exists a maximum coronal temperature 
above which no self consistent solutions exist. 

Since the major physical cause preventing the sandwich con- 
figuration ffom explaining the observed hard spectra is the very 
efficient Compton cooling, we have sought alternative geome- 
tries that do not suffer from a large covering factor. The need for 
a smaller covering factor is also seen ffom the weakness of the 
reflection features in GBHC compared to those in AGN. In our 
previous work we have advocated one such configuration, the 
“sphere+disk” geometry (Dove et al. 1997, see also Gierlirtski 
et al. 1997; Poutanen, Krolik & Ryde 1997). Here, a hot spher- 
ical corona surrounds the black hole, while the cold accretion 
disk which provides the seed photons is situated farther outside. 
This geometry obviously has a smaller covering factor, and in 
fact much higher coronal temperatures are possible (Dove et al. 
1997). Using self consistent numerical models computed for 
the sphere+disk geometry we were able to describe success- 
fully the hard-state spectrum of Cyg X-l measured by RXTE 
over a broad range in energy, 3-200 keV (Dove et al. 1997; pa- 
per I). We derived an optical depth for the spherical corona of 
r = 1 .6 ± 0. 1 and an average temperature of kT = 87 ± 5 keV 
(paper I). 

To facilitate comparison with previous work, and to test 
for the presence of features in the broad-band spectrum, we 
also performed purely phenomenological fits to Cyg X-l. We 
found that the 3 to 200 keV RXTE spectrum could be described 
by an exponentially cut-off power law with a photon index 
T = 145^ q 2 , e-folding energy £* = 162^ keV, plus a devia- 
tion ffom a power law that formally can be modeled as a ther- 
mal blackbody with temperature kT\> b = 1 .2^ ° } keV. Note that 
the above power law index is harder than that typically found, 
which might be an indication that the source was not in its usual 
“hard state” after exiting the “soft state” of 1996 (Cui et al. 
1997a; Zhang et al. 1997 and references therein), which ended 
only a few weeks prior to our observation. No hardening was 
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Table 1: Parameters for the system Cyg X-l/HDE 226868. Values in parentheses are uncertainties in units of the last digit. 



Optical 
Companion 
HDE 226868 

Compact 

Object 

CygX-l 

System 

Parameters 

References 

Position: 





J2000.0: a 



^s^^oo 

1 

<5 



+35° 12'05"82 

1 

galactic: /n 



71.°33 

1 

^ii 



+3?07 

1 

Spectral Type 

09.7 lab 



2 

Teff [K] 

32000 



3 

£b-v 

0.95(7) 



4 

Nh [cm" 2 ] 



6(2) x 10 21 

4,5 

Distance [ kpc] 



2.5(3) 

6 

m v [mag] 

8.84 



7 

B t [mag] 

9.828(22) 



7 

V T [mag] 

9.020(17) 



7 

B-V [mag] 

+0.81 



8 

U-B [mag] 

-0.28 



8 

Luminosity [L 0 ] 

10 54 

10 4 


3,9 

Luminosity [erg/s] 

10 39 

4 - 10 37 



Inclination 



35° 

10 

asini [km] 


5.82(8) x 10 6 


11 

Orbital Period [d] 



5.59974(8) 

11 

Major Axis [R 0 ] 

14.6 

26.3 


10 

vsin / [km/s] 

75.6(10) 



11 

Mass Function [M 0 ] 

0.252(10) 



11 

Mass [M 0 ] 

18 

10 


10 

Radius 

17R© 

30 knit 


3 

Separation [R 0 ] 



41 

10 

Mass loss rate [M 0 /a] 



3x KT 6 

3 

Wind velocity [km/s] 

2100 



10 


, Tiiron et al. (1992), 2 Walbom (1973), 3 Herrero et al. (1995), 4 Wu et al. (1982), 5 Balucirtska-Church et al. (1995), 6 Ninkov, Walker 
& Yang (1987b), 7 ESA (1997), Volume 9, 8 Lutz & Lutz (1972), 9 Liang & Nolan (1984), ,0 see text for further explanation, n Gies & 
Bolton (1982) 

1 Schwarzschild radius 

Table prepared in collaboration with Katja Pottschmidt (Pottschmidt 1997). 
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seen in the HEXTE data from 30 to lOOkeV, indicating that 
any reflection features were formally very weak. [That is, the 
presence of reflection needs to be inferred from more sophis- 
ticated continuum models (Dove et al. 1998; Gierlinski et al. 
1997; Poutanen, Krolik & Ryde 1997), rather than from simple 
exponentially cutoff power laws reflected from cold slabs.] 

2. POWER SPECTRAL DENSITY 
2.1. Techniques and Noise Levels 

The RXTE observations presented here were performed Oc- 
tober, 22, 1996. For the timing analysis we used data obtained 
with the Proportional Counter Array (PCA), the low energy in- 
strument onboard RXTE. The PCA consists of five nearly iden- 
tical Xe proportional counters units (PCU) with a total effective 
area of about 6500cm 1 2 (Jahoda et al. 1996). The data extrac- 
tion was performed using the RXTE standard data analysis soft- 
ware, ftools 4.0. To avoid contamination of the data due to the 
Earth’s X-ray bright limb, only data measured at source eleva- 
tions more than 10° above the spacecraft horizon were used. 

As discussed in paper I, throughout our observations Cyg X- 
1 showed count rates of « 4300 cts s _1 . The total integrated on 
source time for our observation was slightly more than 20 ks. 
For approximately 2 ks two PCU of the PCA were not opera- 
tional, so these data were ignored. We therefore had approxi- 
mately 1 8 ks of high quality timing data 1 . The PCA data modes 
for our observation are listed in Table 2, the formal character- 
istics of the modes can be found in Appendix F to the NASA 
Research Announcement for RXTE (NASA 1997). 

We divided our data into 5 energy bands with 0(700- 
1000 cts s -1 ) each, binned at 2~ 9 s resolution 2 . We also created a 
very high energy band which our energy spectral analysis indi- 
cated should be entirely dominated by background counts. This 
was done to assess the level at which background fluctuations 
could affect our results at lower energies. Thus, we created 
lightcurves for six energy bands in all, labeled 1-6. The prop- 
erties of these energy channels are listed in Table 3. 

We calculated Power Spectral Densities (PSD) in each chan- 
nel by dividing the data into contiguous segments of uniform 
length and time resolution, and performing a Fast Fourier 
Transform (FFT) on each data segment. (Segments with data 
gaps of any duration were ignored.) We performed all the FFTs 
using the reverse discrete transform from EDL (Version 4.0.1), 
which corresponds to 


PSD, 

<|S; I 2 ) = (S)Sj) , (2) 

where the angle brackets indicate an average over data segments 
and , in most cases, over frequency bins. Throughout we have 
used a logarithmic frequency binning wherein we averaged over 
frequencies /— v f+df y with d/// = 0.15. 

We have employed the PSD normalization of Miyamoto 
et al. (1991). In this normalization, one multiplies (\Sj\ 2 ) by 
27l~ 2 T~\ where 71 is the mean count rate (in cts $ -1 ) for the 
energy band of interest. As \Sj\ 2 = |S_,| 2 , we have adopted a 
‘one sided’ normalization. The factor of 2 in the definition of 
the PSD is introduced so that integrating the PSD over positive 
Fourier frequencies and then taking the square root yields the 
normalized root mean square (rms) variability 


rms = 


/ <, 2 )-(,) 2 V /2 
v <*> 2 ) 


(3) 


where here the average is over individual data segments and 
time bins. Note that rms variability is often quoted for the 
PSD integrated over restricted frequency ranges. We shall al- 
ways explicitly state the frequency ranges used in our mea- 
surements. Throughout the rest of the paper we shall take 
Psifj) = 27V 2 T^\S(Jj)\ 2 to signify a PSD as calculated from 
a discrete FFT. When considering two energy channels, P s will 
indicate the PSD of the soft channel, and P h the PSD of the hard 
channel. 

To leading order (cf. Leahy et al. 1983; van der Klis 1989), 
counting noise leads to an additive “white noise” component. 
Its PSD, which we write P*(/), has an amplitude of 2/71, in- 
dependent of Fourier frequency. Its error is 0(1) for a single 
Fourier frequency measured from a single data segment, and is 
distributed as 1 /T2. times x 2 with 2 degrees of freedom 3 * . Instru- 
mental deadtime (cf. Zhang et al. 1995) modifies this slightly. 
By deadtime we mean any process where the detection of an 
event subsequently changes the sensitivity of the detectors [see 
Zhang et al. (1995); Zhang & Jahoda (1996) for a discussion 
of idealized models that can be used to approximate the effects 
seen in RXTE data]. Here we use the deadtime model of Zhang 
et al. (1995), as was successfully applied by Morgan, Remillard 
& Greiner ( 1 997) to RXTE data of GRS 1 9 1 5+ 1 05 . The expec- 
tation value for the PSD of the noise level is different for each 
energy band, and is given by 


AM 

Sj = Y s k txp(2nijk/N) . (1) 

*= o 

Here is the counts in the k? h bin of the lightcurve (consist- 
ing of N evenly spaced time bins), and Sj is the complex dis- 
crete Fourier transform corresponding to the Fourier frequency 
fj - J/N Tj -1 * where T s is the length (in seconds) of each data 
segment and j E [-N/2, N/2]. From 5;, we form the discrete 


P»(f) 




N f - 1 
N f 


7Zp e Td 


+ 7Zp e 7Z v ic 


(?) 


sin(7rr Wt /) 


COS(27Tfj,/) 


(4) 


1 Note that this is slightly less than twice the duration of the observations we used for the spectral analysis of paper 1. Approximately half the standard2f were 
not included in the first data tape that we received; however, all the timing data were included. 

2 We also created lightcurves from one quarter of our 2~ 12 s resolution data to search for high-frequency signatures. No variability in excess of counting noise was 
observed above » 200 Hz. We therefore chose to perform our analyses with 2“* s resolution data. 

3 The fractional uncertainty in the signal component of the PSD is also of 0(1), as measured for an individual Fourier frequency from the FFT of a single data 

segment (cf. Bendat & Piersol 1986; van der Klis 1989). 
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where H e is the total count rate in the energy band of interest, 
1Z pe is the count rate per PCU in the energy band of interest 
(as all five PCU were always on for our selected observations, 
H pe = ^r/5), Uvu is the ‘very large event rate* (mainly due to 
high energy charged particles, and < 1 00 cts s” 1 ), r d is the de- 
tector deadtime (for the PCA we take r d - 10 ps; W. Zhang, Pri- 
vate Communication), r v / f is the deadtime for very large events 
(preset to 70 ps for our observation), Nf is the number of fre- 
quency bins in the FFTs (Nf - N / 2), and tt is the duration of 
the time bins in the input lightcurves. In practice, the terms in- 
volving the cosine and sine above are negligible for the count 
rates and frequencies of interest to us. The above formula was 
compared to the 600-2048 Hz PSDs calculated from our single 
bit data. Assuming that those frequencies were entirely domi- 
nated by counting noise modified by deadtime effects, eq. (4) 
was accurate to better than 2 x \0~ 2 /H e (absolute). Thus eq. (4) 
correctly predicts the deadtime deviation from the uncorrected 
noise level of 2/7Z e to better than 25%. 

We use eq. (4) throughout this work to estimate the deadtime- 
corrected Poisson noise level. All averaged PSD we present 
have had this level subtracted off, as have all quoted rms vari- 
ability levels. The uncertainty in the mean of the PSD is com- 
puted by dividing the uncorrected mean PSD by the square 
root of the total number of independent segments and/or fre- 
quency bins averaged together (cf. Leahy et al. 1983; van der 
Klis 1989). The minimum mean PSD level detectable with our 
observations is approximately the counting noise level from 
eq. (4) divided by the square root of the number of data seg- 
ments and/or frequency bins averaged over. We explicitly show 
this ‘effective noise level’ for all of our average PSDs. The log- 
arithmic binning in frequency results in an approximately f~® 5 
dependence for the effective noise level at high frequency. In 
the next section, we describe our results for Cyg X-l . 

2.2. Results for Cy gnus X~1 

We performed FFTs with three data segment lengths: 1024 s 
(10 data segments), 128s (133 data segments), and 32s (556 


data segments). In Figure 1 we show PSDs calculated from 
these FFTs. As the PSDs for each segmentation of the data 
are statistically correlated, the frequency range 9.7 x 10* 4 - 
1 .3 x 10‘ 2 Hz is computed from the 1024 s data segments, the 
frequency range 7.8 x 10~ 3 -5.9 x 10“ 2 Hz is computed from 
the 128 s data segments, and the frequency range 3.1 x 10~ 2 - 
256 Hz is computed from the 32 s data segments. We use the 
exact same data for calculating the coherence function and time 
lags described below. 

The RXTE deadtime is energy dependent (Zhang & Jahoda 
1996), and in fact the PSD for the background dominated, (45- 
100 keV) band is above the effective noise level. Either eq. (4) 
is inadequate or part of the PSD in this channel may be caused 
by variability of the background or source. However, if the 
10 ps deadtime estimation truly is good to better than ^25% 
at low energy, then our noise-subtracted PSDs are reliable to at 
least 100 Hz. 

For all energy bands, the rms variability in the 10“ 3 -100Hz 
frequency range is « 30%. The PSDs are relatively fiat from 
0.02-0.2 Hz, then fall off at higher frequencies. A portion of 
the higher frequency range has an f~ ] dependence. This be- 
havior is characteristic of the low state of BHC in general (Bel- 
loni & Hasinger 1990a; Belloni & Hasinger 1990b; Miyamoto 
et al. 1992). There is evidence that the PSD steepens again 
at frequencies below 3 x 10“ 3 Hz, especially in the lower en- 
ergy bands. This low frequency noise is partly responsible for 
the greater total rms in the lower energy bands. Low frequency 
noise reported previously has been studied by Pottschmidt et al. 
(1998) and Angelini, White & Stella (1994). There also is a 
weak bump, seen in all energy bands, in at least two frequency 
bins near 0.005 Hz. The rms variability of this region in each 
PSD is « 1 .5%. This feature taken on its own is too weak to be 
called a quasi-periodic oscillation (QPO). A Lomb-Scargle pe- 
riodogram (Lomb 1976; Scargle 1982) of the entire data set also 
shows this feature; however, its significance is only 50-85%. 
(The significance is greatest in the lowest energy band.) As we 
shall discuss in §3 below, however, this region of the PSD, un- 


Table 2: Data modes, with their energy ranges, number of energy channels, and time resolution (always constrained to be a power of 
2 s by the onboard data processors), for our Cyg X-l observation. 


Data Mode 

Energy Range (keV) 

No. Chan. 

Time Res. (s) 

B_2ms_8B_0_35_Q 

0-14.1 

8 

2- y 

SB_250us_0_l 3_2s 

0-5 

1 

2 -l2 

SB_250us_ 1 4_35_2s 

5-14.1 

1 

2-12 

E_125us_64M_36_1s 

14.1-100 

64 

2-13 


Table 3: The energy bands used in this work (labeled 1-6). Also listed are the corresponding PHA channels, their approximate 
energy range, the mean energy of the energy band for our specific observations, and the average count rate in the band. 


Energy Band 

PHA Chan. 

Energy Range (keV) 

Mean Energy (keV) 

Mean Count Rate (cts s 1 ) 

1 

0-10 

0-3.9 

2.5 

780 

2 

11-16 

3 .9-6.0 

4.8 

1040 

3 

17-22 

6.0-8.2 

7.1 

730 

4 

23-35 

8.2-14.1 

10.4 

890 

5 

36-193 

14.1-45 

20.1 

780 

6t 

194-255 

45-100 


70 


t Dominated by background counts. Not used for analysis. Background accounts for « 1% of the counts in band 1, up to te 8% of 
the counts in band 5. Therefore, none of the conclusions in this paper should be strongly influenced by the background. 
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like the surrounding frequency bins, has near unity coherence. 
Low frequency ‘peaks’ in the Cyg X-l PSD previously have 
been observed at frequencies ranging from 0.04-0.07 Hz (An- 
gelini & White 1992; Kouveliotou et al. 1992; Angelini, White 
& Stella 1994; Ubertini et al. 1994; Vikhlinin et al. 1994). For 
the most part these features did not appear as broad peaks with a 
Lorentzian profile. Rather, they appeared as regions in the PSD 
where the power law slope, as a function of Fourier frequency, 
discontinuously changed from positive to negative. Further- 
more, previous observations did not consider the coherence of 
these low frequency features between different energy bands. 

As a phenomenological description of the PSDs, we have 
fit doubly broken power laws over the frequency range 0.02- 
90 Hz. We show an example of one of these fits in Figure 2. 
The reduced \ 2 for these fits ranged from 4 for the lowest en- 
ergy channel to 16 for the highest energy channel. The some- 
what large reduced x 2 values result from the excellent statistics 
of the data, and indicate that doubly broken power laws are only 
qualitative fits to the PSDs. Although the PSD definitely con- 
tain structure not fit by the broken power law, Figure 2 shows 
that they capture the essential features of the 0.02-90 Hz data. 
The overall normalization of the fits decreases slightly with 
increasing energy. The low-frequency power-law index, the 
middle-frequency power-law index, and the two energy breaks 
do not show any significant dependence upon channel energy. 
The high frequency power law index, on the other hand, clearly 
‘hardens’ with increasing energy. This hardening of the PSD 
with increasing energy is evident in Figure 2b. We present all 
of our best fit parameter values in Table 4. 

The results of this doubly broken power law fit are compa- 
rable to the results found by Belloni & Hasinger (1990a) us- 
ing multiple EXOSAT observations. These authors also found 
the PSD to be flat below a break at 0.04-0.4 Hz, then approxi- 
mately oc f~ ] up to a break at 1-6 Hz, and approximately oc f~ 1 
at higher frequencies. There was weak evidence for the high 
frequency power law becoming less steep with increasing en- 
ergy as seen here. Furthermore, the overall rms variability 
ranged from 28%-38%, consistent with the values discussed 
here. Our results and those of Belloni & Hasinger (1990a) are 
also consistent with the trends observed by Ginga , as shown 
in Miyamoto et al. (1992). Belloni & Hasinger (1990a), how- 
ever, did not consider the coherence or time delays between var- 
ious energy bands. As we shall discuss below, these additional 
statistics offer the promise of constraining physical models of 
this system. 
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Fig. 1. — Noise subtracted mean Power Spectral Densities (PSD) for the five energy bands (and one ‘background’ band, see text) 
considered in this work. Top left: (0— 3.9keV). Top right: (3.9— 6.0keV). Middle left: (6.0— 8.2keV). Middle right: (8.2- 14. 1 keV). 
Bottom left: (14.1 -45keV). Bottom right : (45- lOOkeV). For all PSDs, diamonds represent the average PSD with a normalization 
such that the square root of the integrated PSD is the fractional root mean square (rms) variability (see text and Miyamoto et al. 
1991). The solid line represents the “effective noise” level, which is essentially the minimum level of variability that could have 
been detected with our observations. All quoted rms variability levels are for the frequency range 9.7 x 10^*- 100 Hz, except for the 
(45- lOOkeV) band which is the 9.7 x 10" 4 - 10 Hz range. 









8 


RXTE Observation of Cygnus X-l : II. 


Table 4: Parameter values for a doubly broken power law fit to the 0.02 - 90 Hz PSD data. The lowest frequency PSD is fit to a 
function of the form A n f a] . The location of the first break is at 0\ and a 2 is the best fit index of the second power law. The location 
of the second break is at /? 2 and <*3 is the best fit index of the third power law. The reduced Xr values for all the fits are for 44 
degrees of freedom. Uncertainties given are at the 90% level for one interesting parameter (Ax 2 = 2.7). Parameters without errors 
have uncertainties smaller than the last shown decimal place. 


Channel Energy 


c*\ 

A [Hz] 

q 2 

ft [Hz] 


xi 

(0-3.9keV) 
(3.9-6.0keV) 
(6.0-8.2keV) 
(8.2- 14.1keV) 
(14.1 -45keV) 

0.124 ±0.001 
0.116±0.001 
0 1 in-* 005 

u.i ju -qq27 
0 003 

0 0R4^ w* 

0.090S 

0094™ 

0.092$$° 

0-092™ 

A 1 1 A+0 

u liU -fl,Q07 

0.17 

0.17 

0.17±0.01 

0.19±0.01 

0.18 

mm 

2.1 

22$ 

I 9+^:1 
J'rQA 

-1-67$°° 
-l-68$f 
-1.61* 
— 1 48-h503 
* QQ 

-1 44+01* 
1 -^M)00 

4.0 

6.1 
4.4 

9.1 
16.2 




Fig. 2. — Left: The same (0-3.9keV) PSD presented in Fig. 1, except that here we also show the best fit doubly broken power law 
(fit in the range 0.02-90 Hz, reduced x 2 = 4.0). Right: Best fit power law slopes, from doubly broken power law fits, as a function of 
channel energy. Diamonds are the slopes in the w 0.02 -0.2 Hz range, triangles are the slopes in the w 0.2 -2 Hz range, and squares 
are the slopes in the w 2— 90 Hz range. 
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3. COHERENCE FUNCTION 
3.1. The Coherence Function 

Coherence is a lesser known and lesser used statistic than 
the power spectrum. Like the time delay, it is defined for 
two or more concurrent processes (cf. Bendat & Piersol 1986; 
Vaughan & Nowak 1997). It is closely related to the cross cor- 
relation function, and is particularly useful for studying pro- 
cesses whose power density spectra contain features, such as 
QPO or sudden changes in slope, because it provides a means 
of isolating for study a particular Fourier frequency or range of 
frequencies that are strongly correlated across energy bands. 

The coherence function, 7 2 (/), is a Fourier-frequency- 
dependent measure of the degree of linear correlation between 
two concurrent time series. Specifically, it gives the fraction of 
the mean-squared variability at / of one time series that can be 
attributed to, or equivalently predicted from, the other. 

Consider two statistically independent time series, s(t) and 
h(t), measured concurrently. In this work, we are interested in 
cases where s and h are X-ray light curves measured simultane- 
ously in two energy channels. It is (almost) always possible to 
relate s and h by a linear transformation of the form 

hit ) = f t r (t-T)s(T)dr. (5) 

J -co 

The function t r (r) is called the transfer function between s and 
h. In the frequency domain, eq. (5) becomes 

H(f) = Trif)S(f) ( 6 ) 


where £(/), //(/), and T r (f) are the Fourier Transforms of s(t)> 
h(t), and t r (t) respectively. Imagine making a sequence of mea- 
surements of s(t) and h(t), each of finite duration. If T r (J) is the 
same for each measurement of s and /z, the processes are said 
to be perfectly coherent. That is, a measurement of one, say 
.y(f), enables us to predict the other, h(t). The coherence func- 
tion measures the degree to which T r (f) is constant for the data 
segments and Fourier frequencies over which we average. This 
is equivalent to a measure of the degree of linear correlation, in 
the fractional mean-squared sense, between s and h at /. 

The coherence function is defined by 


l<s*0W))l 2 

<is(/)i 2 )<iw(/)iy 


(7) 


where triangular brackets denote an average over an ensemble 
of measurements. The one sigma uncertainty in 7 2 (/), in the 
case of Gaussian statistics, is 


*7 *</)- 


v^7 2 (/)[1~7 2 (/)] 

Yl{f)\y/n 


( 8 ) 


where n is the number of measurements averaged in computing 
7 2 (f) (Bendat & Piersol 1986; Vaughan & Nowak 1997). 

The quantity within triangular brackets in the numerator, 
5* (/)//(/), is the cross power spectrum (sometimes called the 
complex cross spectrum or cross spectral density) between s 
and h. As discussed above, it is a complex quantity whose 
phase is a measure of the shift between the two light curves 
at /. The denominator in eq. (7) is the product of the average 
power spectra. 

An intuitive feel for the coherence function may be gained 
by looking at it geometrically. Let the 7th component of the 


FFT of s (where we now consider an equispaced, discrete-time 
process s*) be Sj , which we write as Sj = The product 

S*Hj is then equal to \S j\\H The coherence function 
is defined as the squared magnitude of the average value of this 
product, | (S* //;)(, normalized by the the product of the indi- 
vidual power spectral densities, (|5y| ){|tf/| }. 





Fig. 3. — a) Vector illustration of the coherence function. Each 
vector is the complex number at a single Fourier fre- 
quency / for a single measurement of A and B. If the two 
processes are coherent, then A*B will always have the same 
phase, as shown in (b). The vectors will be perfectly aligned 
and b® equal t0 If the processes are 

incoherent, then the phase of A*B will be uniformly and ran- 
domly distributed on [ — 7r, 7r], as illustrated in (c). In this case, 
(|A*B| 2 )/((|A| 2 )(|B| 2 )) — ► 0 as the number of measurements 
approaches infinity. 


The process of finding the average value S*Hj may be pic- 
tured as a vector sum in the complex plane. The sum may be 
found geometrically by drawing each term in the sum as a vec- 
tor and adding them head to tail, as in usual vector addition, as 
shown in Figure 3a. If the time series are coherent at frequency 
fj, then the phase difference between them, and hence the di- 
rection of the vector product S*Hj , will be the same for each 
measurement in the sequence, as shown in Figure 3b. Geomet- 
rically, the vectors will all point in the same direction, and the 
sum of the vectors will have the same magnitude as the scalar 
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sum of their magnitudes. In this case, the numerator in eq. (7) 
will equal the denominator, and the coherence is unity. Coher- 
ence can only be computed over an ensemble of measurements. 

In the opposite extreme, the soft and hard channels are com- 
pletely incoherent. The phases of Sj and Hj are unrelated, and 
4>h-4>s will be uniformly and randomly distributed on [-7r,7r]. 
This is just a 2-D random walk, as illustrated in Figure 3c. As 
the number of observations increases, the measured coherence 
goes to zero. 

The light curves s(t ) and h(t) need not look at all alike to be 
perfectly coherent. Indeed, they will look alike only if the time 
difference at each Fourier frequency, Si = (<£>, -<£j)/(27r/), is 
constant, which is equivalent to t r (r) = kS(r-St), where k is a 
scaling factor between s and h. 

An example of near perfect coherence can be found in the 
data studied for this paper. Figure 4 shows lightcurves of Cyg 
X-l in the (0-3.9keV) and (14.1^45keV) energy bands. We 
have divided the data into 128 s segments for measuring co- 
herence. The average direction in the complex plane of the 
Fourier transform of each segment over the Fourier frequency 
range 1.5-2. 5 Hz is shown above the top light curve and be- 
low the bottom, labeled (A) and (B). The product (yTZ?) is 
shown between the light curves. Notice that the phase of the 
Fourier transform varies randomly from segment to segment. 
In contrast, (A*B) remains constant, indicating a fixed phase 
shift between the light curves in the 1. 5-2.5 Hz range. As we 
will show later, the coherence is unity over a broad frequency 
range, even though the phase shift varies by a factor of 4 over 
the same range. 

In practice, measurement noise (in our case Poisson count- 
ing noise) causes the measured coherence to be always smaller 
than unity. Poisson noise contributes a random part to each 
phase measurement, causing the {£(/)* //(/)) to be slightly 
misaligned. For weak signals (relative to the noise), counting 
noise dominates the signal. Poisson noise always dominates 
at sufficiently high frequencies because it has a (nearly) white 
power density spectrum. Coherence is also difficult to measure 
at low frequencies because the number of independent measure- 
ment is smaller. Vaughan & Nowak (1997) discuss methods to 
correct 7 2 and 6~f 2 for counting noise. 

3.2. Coherence in Cyg X-l 

In Figure 5 we present the coherence function for our Cyg 
X-l observation. Specifically, the coherence function and its 
error bars were calculated by applying eq. (8) of Vaughan & 
Nowak (1997) to our data. (Strictly speaking that equation is 
most valid for the high coherence, high signal regime of our 
PSDs, which here is / <> 10 Hz . Improved estimates using 
the more complicated formulae presented in Vaughan & Nowak 
(1997) are unlikely, however, to drastically change our results.) 
All comparisons shown are to the (0-3.9 keV) energy band. As 
discussed by Vaughan & Nowak (1997), the expectation from 
many theoretical models is that the coherence should be less 
than unity. As with previous observations of Cyg X-l (Vaughan 


1991; Vaughan & Nowak 1 997 ; Cui et al . 1 997b), the coherence 
is remarkably close to unity over a wide frequency range, even 
for rather disparate energy bands. Over the range 0.02- 10Hz, 
the coherence is nearly unity. In the two highest energy chan- 
nels, there is a trend for the coherence to drop slightly above 
« 1 Hz. In all energy channels, there is a trend for the coher- 
ence to drop below & 0.02 Hz and above m 10 Hz. For the most 
part, the deviation from unity coherence becomes greater with 
increasing energy. At the lower frequency end, the coherence 
recovers back towards unity near 0.005 Hz, coincident with the 
weak low frequency feature seen in the PSD. 

Further of note is that in the 0.02 -0.2 Hz range, energy band 
5 is actually slightly more coherent than energy band 4 (both as 
compared to energy band 1). Due to the extremely good mea- 
surements in this frequency range, this difference is actually 
significant. The frequency regime 0.02-0. 2Hz corresponds to 
the flat part of the PSDs for energy bands 1-5. 

We highlight these features in Figure 6 by presenting a grey 
scale plot of log[l — 7 2 (/)], which enhances the small differ- 
ences from unity coherence. In this figure, high coherence is 
white, and low coherence is black. The lowest energy chan- 
nel, by definition, is coherent with itself. The general trend to 
lose coherence at both low and high frequencies is evident in 
this plot, as is the trend for this loss to increase with energy. 
The highest energy channel does indeed appear more coherent 
in the 0.02-0.2 Hz range. Note that the high coherence coinci- 
dent with the 0.005 Hz feature is also evident in Figure 6. 



Fig. 4. — Lightcurves (solid lines) for a portion of our Cyg 
X-l data. The top curve is the (14.1— 45keV) energy band, 
shifted by 600 cts s' 1 , while the bottom is for the (0-3.9 keV) 
energy band. Solid vertical lines divide the data into segments 
of 128 s. Arrows represent the directions in the complex plane 
of the FFT and cross power spectral density for each data seg- 
ment, averaged over the frequency range 1.5-2. 5 Hz. Poisson 
noise has not been subtracted; therefore, it is expected that 
(A*)(B) J ( A*B ), as seen for these data. 





Frequency (Hz) 

Fig. 5. — Coherence, as a function of Fourier frequency, 
line is coherence equals unity.) Error bars are calculated from eq. 
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Fig. 6 . — Grey scale representation of the coherence function 
as a function of frequency (vertical axis) and energy band (hor- 
izontal axis). Specifically, there are 256 grey scale levels cor- 
responding to the value of log[l -7 2 if)], where 7 2 (/) is the 
noise-subtracted coherence function for the indicated energy 
band as compared to the (0-3.9keV) energy band. Values of 
7 2 (f) > 0.999 were set to 0.999. White corresponds to the high- 
est coherence levels, black to the lowest coherence levels. 

The deviations to greater than unity for the noise-subtracted 
coherence in Figures 5 and 6 are due to both statistical and 
systematic errors. As the coherence function is a fourth order 
statistic [as compared to a second order statistic like the PSD, 
cf. Bendat & Piersol (1986)], it is very sensitive to uncertain- 
ties in the noise subtraction. With a systematic uncertainty in 
the deadtime as large & 25%, and even with the background 
variability levels as seen in Figure 1, we are still very confi- 
dent of the coherence estimates from « 10 _3 -30Hz, where sys- 
tematic uncertainties are essentially negligible. Between « 3Q- 
100 Hz, our estimates are that systematic uncertainties increase 
the shown error bars by no more than w 30%. The noise- 
subtracted coherence points do show that our noise subtraction 
is not perfect; however, the general trend of coherence loss is 
likely fairly secure to « 100 Hz. 

4. TIME LAGS 

4. 1 . Simple Considerations 

The Fourier time delay, or lag, is also computed for two con- 
current and correlated time series. Like the coherence, it is re- 
lated to the cross correlation, and is a Fourier-frequency de- 
pendent measure of the time delay between the time series (cf. 
Miyamoto & Kitamoto 1989; Miyamoto et al. 1992; van der 
Klis 1989). 

As in §3, let s(t) be a discretely and uniformly sampled “soft 
energy” lightcurve and h(t ) be a discretely and uniformly sam- 
pled “hard energy” lightcurve. These lightcurves, as for the 
PSD (cf. §2) and coherence (§3), are divided into data segments 
of uniform length and an FFT is performed on each segment. 
The Fourier phase lag, <£(/), is the phase of the average cross 
power spectrum. That is, = arg[C(/)], where 

C(/)«<r (f)H(f)) . (9) 

The Fourier time lag is constructed from 4>{f) by dividing 
through by 2irf , i.e. r(f) = <^(/)/2tt/. The Fourier time lag 




(0-3.9 k*V) (3 .9-0.0 k*V) (6 .0-8 .2 k«V) (8.2-14.1 k«V) (14.1-45 k*V) 

Energy Channel 


can be either positive or negative. For our sign convention, a 
positive time lag indicates that the hard light curve lags the soft 
light curve. Note that the Fourier phase is defined on the in- 
terval (-7T, 7r), hence a lag of tt/2 in reality could be a lead of 
37t/2. For a further discussion, see Nowak & Vaughan (1996). 

In general, both the phase lag and time lag are non-constant 
functions of Fourier frequency, /. To get a feel for what this 
can mean, consider the following simple example. Imagine we 
have a source of fluctuations that produces soft X-ray photons 
some distance from a region that also responds to these fluctu- 
ations by producing hard X-ray photons. If the fluctuations can 
propagate from the soft-photon-producing region to the hard- 
photon-producing region, without dispersion, then we expect 
a time delay between the soft and hard photons that is inde- 
pendent of Fourier frequency, /. The time delay at all Fourier 
frequencies will simply be the distance between the soft X-ray 
source and the hard X-ray response divided by the propagation 
speed. This means that the Fourier phase lag, <f>(f), will in- 
crease linearly with / (modulo integer multiples of 2 ; r). 

If the two processes are perfectly correlated (cf. §3), the two 
processes are related by a constant, linear transfer function. The 
time lags observed in Cyg X-l are not independent of Fourier 
frequency; however, the phase lags are roughly independent of 
/. Let us look at a crude but illustrative transfer function that 
produces a constant phase lag. Specifically, we consider a trans- 
fer function, T r (f), with both constant phase and amplitude. 
This simple example is in fact a reasonable approximation to 
Cygnus X-l. In this case, we can write 

H(f) = Aexp(±iA*)S(/) , (10) 

where both A and A <j> are real constants. The plus sign is for 
/ > 0 and the minus sign is for / < 0, because the Fourier trans- 
form of a real time series satisfies H{f) = //*(-/). Taking the 
inverse transform oMexp(±/A<£), the transfer function, t r (r), 
can be written in the time domain as 


r r (r) = A 


r 


(id 


That is, a fraction Acos(A<^>) of the hard variability is exactly 
coincident with the soft variability, while a (typically smaller) 
fraction is delayed from the soft variability with a r~ l “tail” in 
the transfer function. For A ~ 1 and A<f> ~ 0 . 1 radians, most of 
the soft and hard lightcurves are exactly coincident with one an- 
other, and only a small fraction of the hard photons lag behind 
the soft. To illustrate the above points, we present examples of 
both constant (as a function of Fourier frequency) time lags and 
constant phase lags in Figure 7. 

The time delays between hard and soft photons may be due 
to a physical separation between a source and a response, as in 
the first example, or they may be related to the emission pro- 
cess itself. For the case of Comptonization, hard photons un- 
dergo more scatterings than soft photons and therefore are nat- 
urally delayed from soft photons. This latter process occurs on 
very short timescales (on the order of the light crossing time of 
the Comptonizing region) and represents the minimum expected 
time delay between hard and soft X-rays within the context 
of this particular model (Miyamoto & Kitamoto 1989; Miller 
1995; Nowak & Vaughan 1996; Kazanas, Hua & Titarchuk 
1997). On the other hand, the maximum expected time delay 
between hard and soft X-rays will be less than the size of the 
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emitting region divided by the slowest propagation speed in that 
region (whether it be a sound speed or a thermal wave speed, 
etc.). As discussed below, most of the hard and soft photons 
seen in Cyg X-l are nearly temporally coincident, and there- 
fore the (possibly extended) sources of soft and hard X-rays are 
likely to be nearly spatially coincident as well. The more spa- 
tially coincident the soft and hard X-ray sources/responses are, 
the shorter the observed time delays will be. 



Time (arbitrary units) 


Fig. 7. — Gaussian of unit amplitude and variance (solid curve), 
and same Gaussian shifted in time by Ar = 0.5 units (dashed 
line) and in phase by A <j> = 2?r/oAr, (dot-dashed line) where 
/o is the lowest Fourier frequency (here /o = 1 /8). Notice that 
both cause a shift to later times, while the constant phase shift 
also induces an extended tail. 


4.2. Noise Limits on Time Lag Measurements 

A geometrical interpretation of the phase/time lags is help- 
ful in elucidating their nature (cf. §3.1). To understand the 
effects of counting noise on measurements of the time lags, 
we again consider the Fourier components of a measured pro- 
cess as vectors in the complex plane. The vectors consist of 
the sum of a signal part and a noise part. Let the measured 
Fourier components in the soft and hard channels be given 
by S(f) = and H(f) = A h (f)r a, +A^(f). 

Here A s and Ah are complex quantities representing the signal 
that have a constant amplitude, but have phases (i.e. direc- 
tions in the complex plane) that are a function of frequency, 
/. The quantities A^ and A^ are the noise components of 
the Fourier transforms. They have a well defined mean am- 
plitude but a completely random phase. The phases of the 
noise components — and thus their orientation in the complex 
plane — vary randomly from frequency to frequency and from 
data segment to data segment. For simplicity, we shall take 

A,(/)r- =A b (f)r ak = A(J)r a and K*>(/)| = KV)| = 

which is approximately true for our Cyg X-l data. 

As for the coherence function, when we average the CPD 
over discrete samples of the lightcurve and/or over frequency 
bins, we are essentially summing vectors in the complex plane. 
For unity coherence, the signal vectors each have constant di- 
rections set by the phase lags and amplitudes set by the PSD 
(here equal to A 2 f~ 2a ). More generally, if we sum the sig- 
nal over N lightcurve samples and/or frequency bins, the total 


length of the summed signal vectors becomes 

N , ( 12 ) 

where the factor of \J ~f 2 (J) takes account of the coherence 
function (see §3). In our data set, however, the coherence func- 
tion is unity below ~ 10 Hz, and falls off rapidly at higher fre- 
quencies. 

The noise vectors associated with the cross power spectral 
density each have a mean amplitude of 2y/Ps\A\f~ a (i.e. twice 
the product of the amplitude of an individual noise vector with 
the amplitude of an individual signal vector). Their phases, 
however, are completely random. Thus, in the averaging pro- 
cess the sum of these noise vectors looks like a random walk 
in the complex plane (Vaughan 1991, and §3.1 above). The 
root mean square amplitude of the sum of the noise vectors is 
therefore 

s/2NP n \A\r a • (13) 

The fact that this sum has a random phase leads to the uncer- 
tainty in the phase lag determination. In the small angle approx- 
imation (appropriate for large V), the uncertainty of the phase 
lag is then of order the expectation value of the magnitude of 
the sum of the noise vectors divided by the mean magnitude of 
the sum of the signal vectors. We have: 

" fa R ' (l4> 



Fig. 8. — The process of measuring a phase lag in the pres- 
ence and absence of noise, illustrated as a vector sum in the 
complex plane for the case of perfect coherence. Signal vec- 
tors are displayed as solid lines in (a) and (b), and noise vectors 
as dashed lines in (b). For perfect coherence, all signal vec- 
tors are aligned, as seen in (a) and (b). The noise vectors in 
(b) have phases uniformly and randomly distributed on [-?r, tt], 
and contribute component to the vector sum with random phase 
and an amplitude that grows as the square root of the number of 
measurements. In contrast, the amplitude of the sum of signal 
vectors grows as the number of measurements. 

The process of measuring a phase lag in the presence and 
absence of noise is illustrated in Figure 8 for the case of per- 
fect coherence. All signal vectors, displayed with solid lines, 
are aligned. The noise vectors, displayed as dashed lines, have 
phases uniformly and randomly distributed on [-7r,7r]. The 
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noise vectors contribute a random component to the vector sum. 
The noise component has random phase, and an amplitude that 
grows as the square root of the number of measurements. In 
contrast, the amplitude of sum of signal vectors grows as the 
number of measurements. 

Typically, one has total observing time T segmented into 
lightcurves of length T s (yielding a minimum Fourier frequency 
\/T s ). The FFTs are also logarithmically binned over a fre- 
quency interval with fractional width d. The total number of 
Fourier components averaged at frequency / is then just the to- 
tal number of lightcurves, T /T St times the number of frequency 
bins averaged, df /Tf l , yielding N = dTf. For the PSD nor- 
malization used in §2, |/y| « 2/72, independent of frequency. 
We then have 

A m 2 [d7X 7 2 (/)]' 1/2 Up 1 /“' 1/2 • (15) 

For our Cyg X-l data 7 2 (f) a* 1 and a — 0.5 (i.e. PSD a Z” 1 ) 
in the frequency range ~ 0.2-2 Hz; therefore, A <j>(f) is ap- 
proximately constant. Above ~ 2 Hz, o — 0. 7-0.9 and 7 2 (f) 
begins to drop rapidly. It therefore becomes difficult to mea- 
sure phase/time lags above ~ 30 Hz. This situation is not easily 
improved, even with substantially longer integration times. 

As discussed below, the phase lags in Cyg X-l are roughly 
oc /° 3 (i.e. time lags oc 7 ), and A <f> is comparable to the 
noise near 0.1 Hz. As a ^ 0.5, the weak, positive depen- 
dence of phase lag on frequency is what allows us to detect the 
phase/time lags in the 0.1-30 Hz range. Below / ~ 0.1 Hz, 
a ~ 0, and therefore A oc f ~° s , assuming <f> oc /° 3 . The 
phase/time lags therefore become undetectable below & 0. 1 
Hz. Extending our measurements to a decade lower in Fourier 
frequency would require an integration time of order 40 times 
longer. In Figure 9, we present these idealized estimates of the 
noise limit for parameters characteristic of our observations. 



Fig. 9. — The dashed line is the idealized estimate of the 
noise level for time lag measurements typical in this paper [cf. 
eq. (15)]. We have used the double broken power law fit and 
count rate for the (0-3.9 keV) energy band in our estimate, and 
we have taken 7 2 (f) = exp[-(//30 Hz) 2 ]. The solid line is an 
idealized representation of the time lag we see between our low- 
est and highest energy bands. For frequencies below ~ 0. 1 Hz 
and above ~ 30 Hz, lags cannot be measured because of noise 
limitations. 

As discussed by Bendat & Piersol (1986), we can make these 
noise estimates more rigorous. Specifically, the estimate that 
we use for the uncertainty in the phase is 

= N ~^ 2 \ I . 


where N t as before, is the number of lightcurves and/or fre- 
quency bins averaged over and g(f) is defined as 


_ Kg COgCO) I 2 

- <|S(/)| 2 )<|//(/)| 2 ) ’ 


(17) 


In the above, S(f) and H{f) are the Fourier transforms of the 
two light curves to be compared without Poisson noise sub- 
tracted. This formula has the properties discussed above (in- 
cluding taking account of the coherence function, cf. §3), and 
provides a more rigorous estimate of the errors in the phase lag 
than simple error propagation estimates (e.g., Cui et al. 1997b). 
In fact, g(f) is simply the measured coherence of the two light 
curves, without correcting for counting noise. The error in the 
time lag is then simply given by A r(f ) = A<f>(f)/2nf . 


4.3. Time Lags in Cyg X-l 

In Figure 10, we present the time lags vs. Fourier frequency 
for energy bands 2-5 vs. energy band 1 . The measurements 
are above the noise level, as defined by eq. (16), from 0.1-30 
Hz. Within this frequency range, the data are consistent with 
the softest energy band always leading the harder energy bands. 
Those frequency points where the hard leads the soft are domi- 
nated by Poisson noise fluctuations. 

Some general properties of the observations are notable. 
First, the time lags between all channels approximately show 
a power law dependence upon Fourier frequency, r(f) oc /“° 7 , 
although there are significant deviations from a simple power 
law. This is consistent with previous observations of the low 
state of Cyg X-l (Miyamoto & Kitamoto 1989; Miyamoto et al. 
1 992; Crary et al. 1998). As discussed in §4. 1 , such a frequency 
dependence is closer to a constant phase lag (as a function of 
Fourier frequency) than a constant time lag. The time lags pre- 
sented in Fig. 10 correspond to Fourier phases ranging from 
0.01 (low Fourier frequency, energy band 2 vs. energy band 1) 
to 0.2 radians (high Fourier frequency, energy band 5 vs. energy 
band 1). In terms of a constant phase lag transfer function in- 
terpretation [cf. eq. (1 1)], this means that nearly all of the hard 
and soft flux variations occur quasi-simultaneously, with only a 
small fraction of the hard photons lagging the soft photons. 

The f~° 7 dependence of the time lag also implies a large 
dynamic range in delay times, on the order of a decade and a 
half. The longest time lags are 0.05 sec. Measured as a light 
crossing time, this corresponds to a distance of 10 3 GAf/c 2 , for 
M = 10 A/q. Between energy band 5 and energy band 1, the 
shortest measured time lag is « 10~ 3 sec, which corresponds to 
a light crossing radius of 20 GAf/c 2 for Af = 10 A# 0 . It is this 
large dynamic range in characteristic length scales associated 
with the variability, as further discussed in §5, that is difficult 
for theoretical models to address (although see Kazan as, Hua 
& Titarchuk 1997). However, if one were to invoke slower 
“propagation speeds” (see also paper HI), the implied length 
scales would be commensurately shorter. If one were to restrict 
the size scale to 50 GAf /c 2 , the longest time lags would imply 
propagation speeds of 0.05 c. As discussed in paper HI, this is 
uncomfortably slow for both Advection Dominated Accretion 
Flow (ADAF) models, where one expects the matter in the inner 
regions to be nearly in radial free fall (cf. Narayan 1996), and 
Compton corona models, where the sound speed of the corona 
is a sizable fraction of c. 


(16) 
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The time lags generally increase with energy band separa- 
tion. We quantify this statement in §5 below. This increase 
can be reconciled with either the notion that “propagation path” 
from soft source to hard response increases with energy sepa- 
ration, or that the “propagation speed” from soft source to hard 
response decreases with energy separation. We also note that 
the time lags between bands 1 and 5 are consistent with the 
sum of the time lags from bands 1 to 2, bands 2 to 3, bands 3 to 
4, and bands 4 to 5. 


5. PHYSICAL INTERPRETATIONS 
5.1. PSD Shape 

In paper HI, we consider in detail the implications of our tim- 
ing analysis to the specific Comptonization model that we pre- 
sented in paper I. We also consider in more general terms the 
implications of our timing analysis for ADAF models. Specifi- 
cally, we investigate the phenomenology of what we dub “prop- 
agation models”, and derive constraints on likely propagation 
speeds for disturbances in a disk. Here, however, we adopt a 
slightly more general approach. 

We still wish to consider the implications of Comptonization 
models; however, rather than discuss a specific model, we will 
consider those effects that we generally expect to find. As we 
have discussed elsewhere (Nowak & Vaughan 1996, and refer- 


ences therein), if the source of the observed variability is from 
soft photons, perhaps from an accretion disk, that are subse- 
quently upscattered by a Comptonizing medium, then we ex- 
pect a number of effects to be evident. 

First, at high Fourier frequencies, we expect that dispersion 
in the scattering times in the Comptonizing medium will lead 
to an attenuation in the variability amplitude. This attenuation 
will be manifested as a (typically exponential) rollover in the 
PSD on timescales shorter than the typical diffusion timescale 
through the Comptonizing medium (Brainerd & Lamb 1987; 
Kylafis & Klimis 1987; Wijers, van Paradijs & Lewin 1987; 
Stollman et al. 1987; Bussard et al. 1988; Kylafis & Phinney 
1989; Miller & Lamb 1992; Miller 1995; Nowak & Vaughan 
1996). As discussed in Nowak & Vaughan (1996), the sim- 
plest expectation is an exponential rollover in the PSD that is 
oc exp[-(/r) 2 ], where r ~ 2rcy/N\ t% /c. Here, N is the number 
of scatters that the photon has undergone, and is the mean 
free path between scatters. If the optical depth of the Comp- 
tonizing medium is « 1 , then A es « R , where R is the size of the 
medium. 

In the simplest models of Comptonization (cf. Pozdnyakov, 
Sobol & Sunyaev 1983, and references therein), one expects the 
number of scatters N to be of the order 







Fig. 10. — lime lags, as a function of Fourier frequency, for various energy bands vs. the lowest energy band (0-3.9 keV). Diamonds 
represent hard lagging the soft, whereas asterisks represent soft lagging the hard. Error bars are calculated from eq. 16. Measurements 
are above the noise level in roughly the range of 0. 1-30 Hz. 
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where T is the temperature of the Comptonizing medium, E s 
is the energy of the soft input photons, and E 0 is the energy of 
the observed output photons. For a very simple Compton cloud 
model, one therefore expects an exponential cutoff frequency, 
f c = r~ ] , on the order of 


fc w 50 Hz (^OGAZ/c 2 )' 1 (kT/%0 keV) 1/2 

[ln(£„/l keV)-ln(£ J /0.2keV) + 1.6]' /2 ’ 


(19) 


where we have taken M = 10 M§ and where we have chosen 
a characteristic coronal radius and temperature consistent with 
the spectral model of paper I. 

In Figure 1 1 we show the results of fitting an exponentially 
cutoff power law to the PSDs of §2. For all the energy bands, 
we fit functions of the form af b txp[-(f/f c ) 2 ) to the 3-90 Hz 
data. The reduced \ 2 for the fits ranged from 1 .8 for the lowest 
energy band to 5.2 for the highest energy band. Again, reduced 
X 2 as large as these are partly indicative of the extremely small 
error bars associated with the data, as well as the fact that a cut- 
off powerlaw is only a qualitative description of the data. As 
also discussed in §2 and shown in Figure lib, the PSD slope 
“hardens” as one goes to higher energy bands. The best fit 
slopes ranged from -1.74 for the lowest energy band to -1.4 
for the highest energy band. The change in slope is significant 
and appears to be proportional to the log of the band energy. 

Although exponential cutoffs are allowed for all the energy 
bands, fitting pure power laws (as was done in §2) results in 
nearly equally good fits. Based upon eq. (19), we would have 
expected much stronger cutoffs, if the source of the high fre- 
quency variability is intrinsic to the soft input photons 4 . We 
have fit a function of the form of eq. (19) to our best fit cutoff 
frequencies, which we present in Figure 11c. The best fit nu- 
merator was 123 Hz [i.e. (R/5 0 GM/(?Y\kT /80 keV) 1 ! 2 « 
2.5], and the best fit soft input energy was E s = 2 keV. We em- 
phasize that we do not consider this a detection of an exponen- 
tial rollover, but more of an upper limit to the degree of rollover 
(or equivalently a lower limit to the cutoff frequency, f c ). In 
Figure 11c we present the cutoff frequencies that we naively 
would have expected if E s = 0.2 keV (dashed line); however, 
we have kept the numerator as 123 Hz. The expected value of 
the exponential cutoff frequencies are somewhat lower than the 
values permitted by the data. (We further discuss this point in 
paper HI.) 

The simplest explanation (although see Kazanas, Hua & 
Titarchuk 1997 and §5.3 below) for the increase in the PSD 
“hardness” with energy and for the lack of a strong rollover 
in the high energy PSD is that the high frequencies do not 
represent variability intrinsic to the soft input source, but 
rather are the result of direct modulation by fluctuations in the 
Comptonizing medium itself. As we expect the Comptonizing 
medium to be both very hot and concentrated near the compact 
object where the dynamical timescales are the shortest (Cf. pa- 
per I), it is natural to expect the corona to show harder PSD 
slopes at higher energies. 

We hypothesize that the PSD above the break at w2Hz 
is dominated by fluctuations of the corona on dynamical 
timescales. As discussed in Vaughan & Nowak (1997) and in 
§5.4 below, this is also qualitatively consistent with the loss of 
coherence seen at these high frequencies. We further consider 
the characteristic PSD frequencies below. 





Fig. 11. — Top: The same energy band 5 PSD as presented in 
Fig. 1 with an exponentially cutoff power law fit to the 3-90 Hz 
data. Middle: Best fit power slopes to the 3-90 Hz PSD data 
as a function of channel energy. Bottom; Best fit cutoff fre- 
quency (solid line), using eq. (19). The best fit numerator is 
123 Hz with E s = 2 keV. For the dashed line we have fixed the 
numerator of eq. (19) to 123 Hz, but have taken E s = 0.2 keV. 


Note that this statement is true for a relatively spatiaily compact source. In the model of Kazanas, Hua & Titarchuk (1997), the entire PSD is reprocessed white 
noise from a soft input source. The range of power law slopes is achieved by distributing the Compton scatterings over several decades in radius of the system. We 
elaborate upon this point further in §5.3. 
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5.2. Characteristic Timescales of the PSD 

The Compton corona model that we presented in paper I con- 
sisted of a spherical, central hot cloud surrounded by a thin, 
cold disk. TTie temperature at the inner edge of the cold disk 
was taken to be 200 eV, and there was evidence for an iron line 
with equivalent width « 40 eV (cf. § 1 .3). These latter two facts, 
considering that the model was reasonably consistent with the 
data, suggest a radius for the Compton corona of approximately 
50 GM /? if M = 10 M 0 . Substantially smaller radii would im- 
ply a hotter temperature on the inner edge of the disk, which is 
not consistent with the data (Dove et al. 1998). Substantially 
larger radii would imply very weak iron lines, unless a sub- 
stantial fraction of the coronal energy release were also at large 
radii, as in the model of Kazanas, Hua & Titarchuk (1997). 

Taking 50 GM/c 2 as a characteristic radius, we have the fol- 
lowing characteristic timescales for a 10 Mq black hole accret- 
ing at approximately 5% L^dd (cf. Table 1). The light cross- 
ing timescale is of order 2.5 x 10~ 3 sec, which is of order the 
shortest Fourier time lags seen in the system. The dynamical 
timescale (i.e., the Keplerian orbital period) is 0.1 s, or an or- 
bital frequency of 9 Hz. This is of the same order as the location 
of the break in the PSD at « 2 Hz, and is also consistent with 
the Fourier frequencies above which the soft and hard X-ray 
become incoherent with each other. This further supports our 
hypothesis that the variability above ^ 2 Hz is dominated by 
direct modulation of the corona at radii ^ 50 GM/c 2 . We again 
note that there is no evidence for any energy dependence of the 
PSD slopes below 2 Hz. 

The characteristic viscous timescale is of order the dynamical 
timescale divided by the Shakura-Sunyaev a-parameter and by 
the square of the fractional Eddington luminosity (Shakura & 
Sunyaev 1973; Shakura & Sunyaev 1976). Thus this timescale 
is of 0(50-500 s), or equivalently 0(2 x 1 0 -3 — 2 x 1(T 2 Hz), 
for a ~ 0.1-1. Characteristic thermal timescales are of the or- 
der of the dynamical timescale times the a-parameter, i.e. Is 
for a & 0.1. The break in the PSDs at « 0.2 Hz could therefore 
be associated with a characteristic viscous timescale or ther- 
mal timescale. The PSD a f~ ] between 0.2-2 Hz could be in- 
dicative of viscous and/or thermal fluctuations over a range of 
radii & 50 GM/c 2 . Alternatively, one might wish to associate 
this portion of the PSD with dynamical timescale fluctuations 
at radii £ 50 GM/c 2 . 

The feature at « 5 x 10“ 3 Hz is somewhat problematic. This 
frequency is at the lower end of what we would associate with 
a characteristic viscous timescale (unless a is even smaller 
than 0.1), and is too slow for either a thermal or dynami- 
cal timescale unless it is associated with a very large radius 
(R » 3 x 10 4 GM/c 2 for a dynamical timescale). We again 
stress, however, that based upon the PSDs alone, the evidence 
for a low frequency QPO is very weak. 

5.3. Implications of the Time Lags 

The time lags seen in Cyg X-l are often interpreted in terms 
of Comptonization models (cf. Cui et al. 1997b; Kazanas, Hua 
& Titarchuk 1997; Crary et al. 1998 and references therein). 
Taking the simplest Comptonization models, one naturally ex- 
pects the hard photons to lag behind the soft photons as the hard 
photons have undergone more scatters. The expected time de- 
lay is therefore proportional the difference in the average num- 
ber of scatters required to reach the hard energy channel com- 
pared to the soft energy channel times the average distance trav- 
eled between scatters. Again taking the electron scattering op- 


tical depth & 1 so that Rxt A es , the expected time delay is of the 
order 

Ar * (^Mt) ms ' (2o) 

where Eh is the energy of the hard energy channel and E s is the 
energy of the soft energy channel. (As before, we have taken 
M = 10 Mq.) 



Fig. 12. — Time lags as a function of energy compared to the 
lowest energy band, (0-3.9) keV, for three different Fourier fre- 
quencies. Diamonds are for 0.3 Hz, triangles are for 1 Hz, and 
squares are for 10Hz. Lines are the best fit functions of the 
form of eq. (20). 

In Figure 12 we present the time delay, compared to the low- 
est energy band, as a function of energy for three frequencies, 
0.3 Hz, 1 Hz, and 10 Hz. As has been noted recently for other 
observations (Cui et al. 1997b), the time delays do indeed in- 
crease nearly linearly with the log of the band energy. We have 
fit a function of the form of eq. (20) to these data, and find the 
best fit intercept to be given for Es « 3.0-3.2keV, in reasonable 
agreement with the average energy of 2.5 keV for the lowest en- 
ergy channel. (We did not include a 0 sec time delay between 
the lowest energy band and itself in the fits, and the y-intercept 
was a free parameter.) Other authors (e.g. Cui et al. 1997b; 
Crary et al. 1998) recently have claimed that this previously ob- 
served logarithmic energy dependence of the time lag is strong 
evidence for the presence of a Compton corona; however, rarely 
do they note the fact that although the functional form agrees, 
the best fit slopes are apparently functions of Fourier frequency. 

Specifically, if we set JcT = 80 keV in eq. (20), then we find 
best fit values of R = 230 GM/c 2 for / = 0.3 Hz, R = 85 GM/c 2 
for /= 1 Hz, and = 17 GM/c 2 for / = 10 Hz. Any individ- 
ual Fourier frequency agrees well with the simple expectations 
of Comptonization theory; however, taken as a whole the time 
delays do not agree with the simplest Comptonization mod- 
els. This fact was originally noted by Miyamoto & Kitamoto 
(1989). Two possibilities have been suggested for maintain- 
ing the logarithmic energy dependence, yet creating the large 
dynamical range in time delays, within the context of Comp- 
tonization. The first is the possibility that the time delays are 
inherent to the source of soft input photons and are merely “re- 
processed” by the Comptonizing medium (Miller 1995; Nowak 
& Vaughan 1996). We explore this possibility in detail in paper 

m. 
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The second suggestion is that the dynamic range in time de- 
lays reflects a large range in radii of a Comptonizing medium 
(Kazanas, Hua & Titarchuk 1997). Specifically, Kazanas, Hua 
& Titarchuk (1997) posit a uniform temperature corona with 
essentially equal optical depth per decade of radius (i.e., elec- 
tron density approximately oc ) that extends for roughly three 
decades in radius. A white noise (i.e., flat PSD) input spectrum 
is reprocessed into a power law PSD that shows time delays 
comparable to those shown in Figure 12. Specifically, the time 
delays, at a given Fourier frequency, show a logarithmic energy 
dependence, and higher Fourier frequencies show shorter time 
delays. This latter attribute is due to the fact that any portion of 
the lightcurve that comes from photons that have scattered over 
large distances has its high frequency variability strongly atten- 
uated [cf. eq. (19)]. Any observed high frequency variability is 
only sampling those photons that have scattered over short dis- 
tances, and therefore it will show relatively short time delays 
between soft and hard photons. Furthermore, the r“ ! density 
profile means that there is roughly equal probability of scatter- 
ing on length scales in any given decade of radius, which leads 
to the observed dynamic range in the time delays. Although the 
model of Kazanas, Hua & Titarchuk (1997) is consistent with 
the data, it is difficult to understand how such a configuration 
would evolve. It implies that the integrated electron energy goes 
as r 2 , i.e. most of the electron thermal energy is at the largest 
radii. For the Compton corona model we discuss in paper I, the 
bulk of the accretion and electron thermal energy is expected to 
be at radii ;$ 50 GM/c 2 . 

5.4. Implications of the Coherence 

The physical implications of the coherence function are dis- 
cussed extensively in Vaughan & Nowak (1997). Here we point 
out some of the specific implications for Comptonization mod- 
els. If the low frequency variability is reprocessed from the soft 
input photons, the near unity coherence between « 0.02- 10 Hz 
indicates that on these timescales the Comptonizing medium 
must appear as a static medium. [As discussed by Nowak & 
Vaughan (1996), a static corona is a linear reprocessor.] If, on 
the other hand, the corona directly creates the low frequency 
variability, then it must do so in a global, linear fashion, despite 
the fairly large observed variability amplitudes. 

As discussed by Vaughan & Nowak (1997), the loss of high 
frequency coherence is consistent with the presence of mul- 
tiple, independent sources and responses. As the frequen- 
cies associated with the dynamical timescales are £ 9 Hz for 
R £ 50 GAf/c 2 and M = 10 A/ 0 , the high frequency coherence 
loss and the “hardening” of the high frequency PSD with energy 
is consistent with the hypothesis that we are directly viewing 
the formation timescales and processes of the corona. 

The loss of coherence below <> 0.02 Hz is consistent with 
viscous timescales for R ^ 50 GAf/c 2 . Perhaps on these low 
frequency timescales we are seeing the disk ‘global’ structure 
becoming disassociated from the corona ‘global’ structure. 

We again note that although there seems to be a general trend 
for the coherence to drop for frequencies £> 0.02 Hz, there ap- 
parently is a recovery at 0.005 Hz, coincident with the weak 
feature in the PSD. In general, we expect any global oscillation 
to be fairly coherent. The coherence associated with this feature 
makes it appear to be more significant than implied by the rather 
low 50-85% value suggested by a Lomb-Scargle periodogram 
(cf. §2). Unfortunately, we are not aware of any direct way 
of assessing the statistical significance of this coherence peak 


coincident with the weak PSD feature. 

6. SUMMARY 

We have presented timing analysis for a 20 ks RXTE obser- 
vation of Cygnus X-l, taken approximately two weeks after 
Cyg X-l returned to its hard/low state from its soft/high state. It 
is possible that some of the properties of Cyg X-l were slightly 
different from their “standard” hard/low state values (e.g. the 
10-200 keV photon index). However, many of the timing prop- 
erties were similar to previous observations (Miyamoto & Ki- 
tamoto 1989; Belloni & Hasinger 1990a; Belloni & Hasinger 
1990b; Miyamoto et al. 1992). 

We have presented detailed discussions of all the timing anal- 
yses used in this work, including discussions of errors and noise 
limits. We are very confident of our results for the PSD from 
10 _3 -100Hz, for the coherence from 10“ 3 -30Hz, and for the 
time delay from 0. 1-30 Hz. The coherence function is the most 
sensitive to systematic errors (i.e., deadtime); however, the gen- 
eral trend of coherence loss out to « 100 Hz is secure even with 
uncertainties as large as 25% in the deadtime. ( Ginga system- 
atic uncertainties were such that one could not be confident in 
coherence trends above « 1 Hz; cf. Vaughan 1991; Vaughan & 
Nowak 1997.) 

Some of the notable features of our observations were the 
0.005 Hz feature and the high frequency ‘hardening’ with en- 
ergy seen in the PSD. The former, if real, could be indicative of 
an oscillation on viscous timescales, whereas the latter may be 
indicative of dynamical timescale variability produced directly 
by the corona. The coherence function lends support to both 
of these interpretations. For the former, the coherence is seen 
to recover towards unity, whereas for the latter, the coherence 
is seen to drop in accord with theoretical expectations of mul- 
tiple flaring regions (Vaughan & Nowak 1997). The coherence 
was also seen to be nearly unity from 0.02-10Hz. Curiously, 
the 0.01-0.2 Hz coherence between the lowest and highest en- 
ergy bands was seen to be slightly greater than the coherence 
between the lowest and second highest energy bands. 

We also made comparisons of our timing analysis to the ex- 
pectations of simple Comptonization models. The high fre- 
quency PSD in the highest energy channels was not as atten- 
uated as expected for simple corona models, again lending cre- 
dence to the hypothesis that we are seeing direct modulation 
by a fluctuating corona. Although the time lags show the ex- 
pected logarithmic dependence upon energy, their increase with 
decreasing Fourier frequency is contrary to the simplest expec- 
tations for a relatively compact, uniform corona. We explore 
this point further in a companion paper (paper EQ). 

Both our energy spectral analysis of paper I and the timing 
analysis presented here show that the simplest models are in- 
adequate for describing the high quality RXTE data. Due to 
its broad spectral coverage, RXTE allows us to compare and 
discriminate among detailed, self-consistent models for the en- 
ergy spectrum. Coupling these detailed models (as we do in pa- 
per HI) with RXTE’s superior timing capabilities places further 
constraints on their structure, and more importantly, on their 
dynamics. Only those models in which both the spectral and 
temporal data can be explained should be considered valid can- 
didates for Cygnus X-l . 
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